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Abstract 

The gravity-driven flow along an annular topological defect (string) with transversely 
corrugations is investigated by using the verified transition-rate model and boundary pertur- 
bation method. We found that for certain activation volume and energy there exists possible 
frictionless states which might be associated with the missing momentum of inertia or dark 
matter. 
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1 Introduction 

Quite recently Vitelli et al. suggested that topological defects in the cosmic shear can be used as 
a probe of the gravitational potential generated by the lensing mass fluctuations on large scales 
[1] based on the facts that shear fields due to weak gravitational lensing have characteristic co- 
herent patterns. They described the topological defects in shear fields in terms of the curvature 
of the surface described by the lensing potential. In that paper they explored the connection 
between the theory of topological defects and the spatial patterns of shear fields due to weak 
gravitational lensing. 

The starting point of their approach rests on an analogy between gravitational lensing shear 
fields, as a probe of structure formation on cosmological scales, and the anisotropic optical or 
mechanical response of materials, as a probe of their inhomogeneous structure on microscopic 
scales. As an illustration, the topological defects in the local shear field of an elastic medium 
reflect the external deformations applied to the solid. Similarly for thin liquid crystal films con- 
fined on a curved substrate, the density of topological defects depends on the inhomogeneous 
curvature of the underlying surface [1] . This may also allow them to infer how the dark matter 
is concentrated around galaxies and galaxy clusters, as well as providing a testing ground for 
dark energy and modified gravity theories [2,3]. 

Above mentioned or borrowed analogy is one motivation for our present study. The other mo- 
tivation is related to the possible dark matter associated with possible superfluidity formation 
after shear-thinning. While superflow in a state of matter possessing a shear modulus might 
initially seem untenable, experimental claims for precisely this phenomenon in solid 4 He now 
abound [4]. Reported in the experiments of Kim and Chan [5] was a dramatic change below 200 
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mK in the period of a torsional oscillator containing solid 4 He. Because superfluids come out of 
equilibrium and detach from the walls of the rotated container, they are expected to give rise 
to a period shift in such a geometry, assuming, of course, the rotation velocity is less than the 
critical velocity to create a vortex. The result is a missing moment of inertia (MMI) [4] and 
hence the period of oscillation decreases. The magnitude of the MMI is a direct measure of the 
superfluid fraction. However, the present author likes to link this MMI which occurs as there is 
formation of superfluidity with the possible formation of dark matter. 

Meanwhile, researchers have been interested in the question of how matter responds to an 
external mechanical load. External loads cause transport, in Newtonian or various types of 
non-Newtonian ways. Amorphous matter, composed of polymers, metals, or ceramics, can de- 
form under mechanical loads, and the nature of the response to loads often dictates the choice 
of matter in various applications. The nature of all of these responses depends on both the 
temperature and loading rate. 

To the best knowledge of the author, the simplest model that makes a prediction for the rate 
and temperature dependence of shear yielding is the rate-state model of stress-biased thermal 
activation [6-8]. Structural rearrangement is associated with a single energy barrier E that is 
lowered or raised linearly by an applied stress a : R± = uq exp[— E/(kBT)] exp[±aV* /(hsT)], 
where ks is the Boltzmann constant, v$ is an attempt frequency and V* is a constant called 
the 'activation volume'. In amorphous matter, the transition rates are negligible at zero stress. 
Thus, at finite stress one needs to consider only the rate R + of transitions in the direction aided 
by stress. 

The linear dependence will always correctly describe small changes in the barrier height, since 
it is simply the first term in the Taylor expansion of the barrier height as a function of load. 
It is thus appropriate when the barrier height changes only slightly before the system escapes 
the local energy minimum. This situation occurs at higher temperatures; for example, New- 
tonian transport is obtained in the rate-state model in the limit where the system experiences 
only small changes in the barrier height before thermally escaping the energy minimum. As the 
temperature decreases, larger changes in the barrier height occur before the system escapes the 
energy minimum (giving rise to, for example, non-Newtonian transport). In this regime, the 
linear dependence is not necessarily appropriate, and can lead to inaccurate modeling. To be 
precise, at low shear rates (7 < j c ), the system behaves as a power law shear-thinning material 
while, at high shear rates, the stress varies affinely with the shear rate. These two regimes 
correspond to two stable branches of stationary states, for which data obtained by imposing 
either a or 7 exactly superpose. 

In this short paper, motivated by the analogy used in [1], we shall adopt the verified transition- 
rate-state model [6-8] to study the gravity-driven transport of cosmic textures (presumed to be 
amorphous) within a corrugated annular (cosmic) string. The possible nearly frictionless states 
due to strong shear-thinning will be relevant to the dark matter formation as mentioned above 
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(considering the MMI [4]). To obtain the law of shear-thinning matter for explaining the too 
rapid annealing at the earliest time, because the relaxation at the beginning was steeper than 
could be explained by the bimolecular law, a hyperbolic sine law between the shear (strain) 
rate : 7 and shear stress : r was proposed and the close agreement with experimental data was 
obtained. This model has sound physical foundation from the thermal activation process [6-8] (a 
kind of (quantum) tunneling which relates to the matter rearranging by surmounting a potential 
energy barrier was discussed therein). With this model we can associate the (shear-thinning) 
fluid with the momentum transfer between neighboring atomic clusters on the microscopic scale 
and reveals the atomic interaction in the relaxation of flow with dissipation (the momentum 
transfer depends on the activation (shear) volume : V* = Vh which is associated with the center 
distance between atoms and is equal to fceT/ro (T is temperature in Kelvin, and to a constant 
with the dimension of stress). 

To consider the more realistic but complicated boundary conditions in the walls of the annular 
(cosmic) string, however, we will use the boundary perturbation technique [10] to handle the 
presumed wavy-roughness along the walls of the annular (cosmic) string. To obtain the analyti- 
cal and approximate solutions, here, the roughness is only introduced in the radial or transverse 
direction. The relevant boundary conditions along the wavy-rough surfaces will be prescribed 
below. We shall describe our approach after this section : Introduction with the focus upon the 
boundary perturbation method. The approximate expression of the transport is then demon- 
strated at the end. Finally, we will illustrate our results into two figures and give discussions 
therein. 

2 Theoretical Formulations 

We shall consider a steady transport of the (shear-thinning) amorphous matter in a wavy-rough 
annular (cosmic) string of r\ (mean-averaged inner radius) with the inner wall being a fixed 
wavy-rough surface : r = r\ + e sin(fc^ + (3) and T2 (mean- averaged outer radius) with the outer 
wall being a fixed wavy-rough surface : r = r2 + e sin(/c#), where e is the amplitude of the (wavy) 
roughness, (3 is the phase shift between two walls, and the roughness wave number : k = 2ir/L 
(L is the wavelength of the surface modulation in transverse direction). 

Firstly, this amorphous matter (composed of cosmic textures) can be expressed as [6-8] 7 = 
70 sinh(r/ro), where 7 is the shear rate, r is the shear stress, and 7o(= C^ksT exp(—AE/kBT)/h) 
is with the dimension of the shear rate; here = 2VyjV m is a constant relating rate of strain to 
the jump frequency (V/, = A2A3A, V m = A2A3A1, A2A3 is the cross-section of the transport unit 
on which the shear stress acts, A is the distance jumped on each relaxation, Ai is the perpen- 
dicular distance between two neighboring layers of particles sliding past each other), accounting 
for the interchain co-operation required, h is the Planck constant, AE is the activation energy. 
In fact, the force balance gives the shear stress at a radius r as r = —{r5Q)/2 [6-8]. 5Q is the 
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net effective gravity forcing along the transport (or tube-axis : z-axis) direction (considering dz 
element) . 

Introducing the forcing parameter $ = — (r2/2To)5Q then we have 7 = 70 sinh(<I>r/r2). As 
7 = —du/dr (u is the velocity of the transport in the longitudinal (z-)direction of the annular 
(cosmic) string), after integration, we obtain 

u = u s + ^p^[cosh<& — cosh(— )], (1) 

here, u s (= u s n p ) is the velocity over the (inner or outer) surface of the annular (cosmic) string, 
which is determined by the boundary condition. We noticed that a general boundary condition 
for transport over a solid surface [9] was 

Su = L° 7 (l - |)~ 1/2 , (2) 

where 5u is the velocity jump over the solid surface, L° s is a constant slip length, j c is the critical 
shear rate at which the slip length diverges. The slip (velocity) boundary condition above 
(related to the slip length) is closely linked to the mean free path of the particles together with a 
geometry-dependent factor (it is the quantum-mechanical scattering of Bogoliubov quasiparticles 
which is responsible for the loss of transverse momentum transfer to the container walls [10]). 
The value of 7 C is a function of the corrugation of interfacial energy. 

With the slip boundary condition [9], we can derive the velocity fields and transport rates 
along the wavy-rough annular (cosmic) string below using the verified boundary perturbation 
technique [11] and dimensionless analysis. We firstly select to be the characteristic length 
scale and set r' = r / LP S , R\ = r±/L®, R2 = r^/L^, e' = e/L^. After this, for simplicity, we drop 
all the primes. It means, now, r, R±, R2 and e become dimensionless (<3? and 7 also follow). 
The wavy boundaries are prescribed as r = R2 + esm(k6) and r = R± + esm(k9 + (3) and the 
presumed steady transport is along the z-direction (microannulus-axis direction). 

2.1 Boundary Perturbation 

Along the outer boundary (the same treatment below could also be applied to the inner bound- 
ary), we have 7 = (du)/(dn)\ on wa \\ s - Here, n means the normal. Let u be expanded in e 

u = u + eui + e 2 u 2 H , 

and on the boundary, we expand u(ro + edr, 6(= 6q)) into 

dv^ 1 

u(r,6)\( ro+edr>eo ) = u(r ,9) + e[dr u r (r ,9)} + e 2 [—u rr (r ,6)] H = 

7i?2 

K/ ip + -|-[cosh$-cosh(— )]} | on walls , r = R 1 ,R 2 ; (3) 



where 



u slip\on walls - ^siM 1 - -r-) 1/2 ]} I on walls' ( 4 ) 

7c 
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Now, on the outer wall (cf. [11]) 

du _ V(r — R 2 — esmikO)) . 9 k 2 ,„„,, i. 

cos(/c6»)u e j (jR2+£(irie )] = u 0r |ii 2 +e[ui r |fl 2 + u 0rr \ R2 sin(fc0)- 
k lk 2 

—uq 6 \ R2 cos(k6)} +e 2 [--— cos 2 (k6)u 0r \ R2 + u 2t \r 2 +u lrr \ R2 sm(k9)+ 



k 

2 -utrr V-- / 



TT^^Iiia sin 2 (A;6') - 3 cos(/c6')(niJ jR2 + u 0gr \r 2 sin(fc0))] +0(e 3 ). (5) 



Considering LP S ~ Ri,R2 3> e case, we also presume sinh<i> <C j c /jo- With equations (1) and 
(5), using the definition of 7, we can derive the velocity field (u) up to the second order : 

u(r,9) = -( J R 2 7o/^){cosh($r/i? 2 ) - cosh$ [1 + e 2 $ 2 sin 2 (A;0)/(2^)] + 

e$sinh$ sin(fc0)/i? 2 } + u s / i p| r= j ?2+esin ( fce ). 

The key point is to firstly obtain the slip velocity along the boundaries or surfaces. After lengthy 
mathematical manipulations, we obtain the velocity fields (up to the second order) and then we 
can integrate them with respect to the cross-section to get the transport (volume flow) rate (Q, 
also up to the second order here) : 

Q= / u(r, 9)rdrd9 = Qq + e Q Po + e 2 Q P2 . 

JO JRi+esin(A:6»+/3) 

In fact, the approximate (up to the second order) net transport (volume flow) rate reads : 

Q = 7Tj {L° s (R 2 2 - R\) sinh $ (1 - $^)-V2 + ^[(#2 _ CQsh $ _ h R 2 sinh $ _ 

7c/ 7o 3> <P 

i^sim!^^)) + (cosh$ - cosh($|i))]} + e 2 {|u sHpo (i? 2 - i? 2 )+ 

L^ 70 sinh$(l + *^){-k 2 + <1> 2 )[1 - A 2 ] + ^70^1 sinh(|i$) - R 2 smh<f>}- 

4 7c/7o -n-2 ^ 

7T i? 2 r , T i?i 7T 3-tt-. fli, 

-70 — [cosh $ -cosh($— )] + -70$ cosh <I>[i? 2 - — -]+ 
2 9 R 2 4 R 2 

7T7o{[sinh<l> + L°cosh<I>(l + S ^^ )](R 2 - Ri cos /?)} + —70-^ cosh$+ 

7c/7o 2 3> 

^n^ T 9 cosh<D r ,JJi,, n , T 



Here, 



^ Po =^7o[sinh$ ( l-^)-V2]. (?) 
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3 Results and Discussions 

We firstly check the roughness effect (or combination of curvature and confinement effects [12- 
13]) upon the gravity-driven transport via strongly shearing because there are no available 
experimental data and numerical simulations for the same geometric configuration (annular 
(cosmic) string with wavy corrugations in transverse direction). With a series of forcings (due 
to imposed gravity forcings) : $ = R.2(5G) / (2tq) , we can determine the enhanced shear rates 
(dj/dt) due to gravity forcings. From equation (5), we have (up to the first order) 



The parameters are fixed below (the orientation effect : sin(A;#) is fixed here). r2 (the mean 



can be tuned easily. The effect of wavy-roughness is significant once the forcing (<&) is rather 
large (the maximum is of the order of magnitude of e[<I> tanh(<E>)/i?2])- 

If we select a (fixed) temperature, then from the expression of To, we can obtain the shear stress 
r corresponding to above gravity forcings (<£) : 



There is no doubt that the orientation effect (9) is also present for the amorphous matter. For 
illustration below, we only consider the maximum case : |sin(fc#)| = 1. We shall demonstrate 
our transport results below. The wave number of roughness in transverse direction is fixed to 
be 10 (presumed to be the same for both walls of the annular (cosmic) string) here. 
Now, we start to examine the temperature effect. We fix the forcing $ to be 1 as its effect is of 
the order O(l) for the shear rate (cf. Fig. 2). As the gravity forcing (5Q) might depend on the 
temperature (5Q = 2to&/R2, tq = tq(T), V*(= Vh) is presumed to be temperature independent 
here for simplicity). Note that, according to [7], V* = 3V5~f/2 for certain matter during an 
activation event [7], where V is the deformation volume, £7 is the increment of shear strain. 
As the primary interest of present study is related to the possible phase transition [14-16] or 
formation of superfluidity (presumed to be relevant to the formation of dark matter mentioned 
in Introduction) due to strong shearing, we shall present our main results in the following. We 
performed intensive calculations or manipulations of related physical and geometric parameters, 
considering a hot big-bang universe [14] and examine what happens as it expands and cools 
through the transition temperature T c . The selected temperature range and the activation an- 

ergy follows this reasoning. Note that in unified models of weak and electromagnetic interactions 

1/2 

T c is of the order of the square root of the Fermi coupling constant [14], Gp , i.e. a few hundred 
GeV. Thus the transition occurs when the universe is aged between and 10 -12 seconds and 
far above nuclear densities [14,17]. One possible superfluidity formation regime is demonstrated 
in Fig. 1. The activation energy (AE) is 10~ 10 Joule (~ 1 GeV) and the activation volume is 
10~ 9 m 3 . In fact, all the results shown in this figure depend on 70 and are thus very sensitive 




(8) 



outer radius) is selected as the same as the slip length L° s . The amplitude of wavy roughness 



r = r sinh [sinh($) + e sin(/c<9) — cosh($)] . 

R2 



(9) 
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to AE. Here = 2 and the sudden jump of the shear stress (directly linked to the friction) 
occurring around T ~ 10 11 °K could be the transition temperature for the selected AE and 
There is a sudden friction drop around two orders of magnitude below T ~ 10 11 °K and it is 
almost frictionless below T ~ 10 10 °K. If we borrow the analogy from the MMI [4] (the rotational 
properties of a superfluid as well as a supersolid in which some of the particles remain still while 
the rest of them rotate with the container) we can identify the formation of superfluidity as the 
possible dark matter formation as the mass is missing within this regime. 

The possible reasoning for this formation can be illustrated in Fig. 2. It could be due to the 
strong shearing driven by larger gravity forcings along a confined (cosmic) string or topological 
defect. The shear-thinning (the viscosity diminishes with increasing shear rate) reduces the vis- 
cosity significantly. One possible outcome for almost vanishing viscosity is the nearly frictionless 
transport. It seems to us that the formation of dark matter is a genuinely dynamic effect. As to 
the question whether the internal structure of defects lead to persistent currents in their cores 
[17] is yet open. We shall investigate the evolution of the (cosmic) string [14,17] in the future. 

References 

[1] V. Vitelli, B. Jain, and R.D. Kamien, Topological defects in gravitational lensing shear 
fields, J. Cosmology Astroparticle Phys. (2009) JCAP09(2009)034. 

[2] Y. Mellier, Probing the Universe with Weak Lensing, Annu. Rev. Astron. Astrophys. 37 
(1999) 127. 

[3] H. Hoekstra and B. Jain, Weak Gravitational Lensing and Its Cosmological Applications, 
Annual Review of Nuclear and Particle Systems 58 (2008) 99. 

[4] P. Phillips and A.V. Balatsky, Cracking the superdolid, Science 316 (2007) 1435. 

[5] E. Kim and M.H.W. Chan, Probable observation of a supersolid helium phase, Nature 427 
(2004) 225. E. Kim and M.H.W. Chan, Observation of Superflow in Solid Helium, Science 
305 (2004) 1941. 

[6] R.E. Powell and H. Eyring, Mechanism for the relaxation theory of viscosity, Nature 154 
(1944) 427. 

[7] S. Glasstone, K.J. Laidler and H. Eyring, The Theory of Rate Processes (McGraw-Hill, New 
York, 1941). 

[8] H. Eyring, D. Henderson, B.J. Stover and E.M. Eyring, Statistical Mechanics and Dynamics 
(Wiley, New York, 1964). 

[9] Z.K.-H. Chu, Preprint (2007). 



8 



[10] D. Einzel and J.M. Parpia, Slip of quantum fluids, J. Low Temp. Phys. 109 (1997) 1. 

[11] W. K.-H. Chu, Stokes slip flow between corrugated walls, ZAMP 47 (1996) 591. 

[12] V. Vitelli and A.M. Turner, Anomalous coupling between topological defects and curvature, 
(2004) |cond-mat/0406329| . 

[13] P. Gaete and C. Wotzasek, Duality and confinement in D=3 models driven by condensation 
of topological defects, [hep-th / 0503070 (2005). P. Gaete and C. Wotzasek, On condensation 
of topological defects and confinement, (2004) hep-th/040 8115| . 

[14] T.W.B. Kibble, Topology of cosmic domains and strings, J. Phys. A : Math. General 9 
(1976) 1387. 

[15] B. Damski, The simplest quantum model supporting the Kibble-Zurek mechanism of 
topological defect production: Landau-Zener transitions from a new perspective, (2004) 
[cond-mat/0411004| . 

[16] M. Uhlmann, R. Schiitzhold and U.R. Fischer, O(N) symetry-breaking quantum quench - 
topological defects versus quasiparticles, (2009) [arXiv:0905.0877]. 

[17] A. Gangui, Cosmology from Topological Defects, (2003) |astro-ph/ 0303504 1 . 



9 



10° 



10° - 



U3 

in 
<u 

L 

J. 10 

TO 

<U 
.c 
CO 



10" 



10" 



❖ o 



- Nearly Frictionless 



10' 



o o 



o o 



Activation energy =10" Joule 



*=1 



io 11 io 12 io 13 

Temperature ( K) 



10" 



Fig. 1. Comparison of calculated (shear) stresses using an activation energy 10~ 10 J or 
~ 1 GeV. There is a sharp decrease of shear stress around T ~ 10 11 K. Below 10 10 K, 
the transport of amorphous matter is nearly frictionless. 
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Fig. 2. Increasing shear causes a local energy minimum to flatten until it disappears 
(energy barrier removal or quantum-like tunneling). The structural contribution 
to the shear stress is referred to shear-thinning. 



